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摘 要 : 冰川 反照 率 作为 影响 辆 射 收 支 最 为 敏感 的 要 素 , 其 波动 变化 将 直接 影响 冰川 能 量 一 物质 平衡 的 变化 。 利 


用 遥感 手段 研究 冰川 反照 率 特征 和 时 空 分 布 规律 ,对 于 估算 冰川 消融 乃至 下 游 河 川 径流 具有 重要 参考 依据 。 本 文 
多 MOD10A1 逐 日 反照 率 产品 和 Landsat TM/ETM+ 影 像 反 照 率 反 演 结果 为 主要 数据 源 , 用 实测 数据 进行 Landsat TM/ 


ETM+ 反 演 反照 率 和 MODI10A1 产 品 反照 率 的 误差 分 析 ; 进 而 分 析 黑 河流 域 典 型 冰川 2000—2018 4E MOD10A1 产品 


反照 率 的 时 空 变化 特征 。 


结果 表明 ,MOD10A1 与 自动 气象 站 测 得 的 实测 反照 率 差 异 较 大 ,但 两 者 随时 间 的 变化 趋 
势 一 致 ;MOD10A1 5j Landsat TMETM+ 反 演 结果 差异 较 小 ,在 黑河 流域 冰川 反照 率 的 研究 中 具有 较 好 的 适用 性 。 


在 年 际 尺度 上 ,4 条 冰川 反照 率 均 呈 现 微 弱 的 下 降 趋 势 , 空 间 尺 度 上 , 八 一 冰川 年 平均 反照 率 最 高 , 羊 龙 河 五 号 冰川 
次 之 ,七 一 冰川 次 低 , 十 一 冰川 最 低 ; 消 融 季 冰川 反照 率 羊 龙 河 五 号 冰川 最 高 ,七 一 冰川 次 之 , 八 一 冰川 次 低 , 十 一 


冰川 最 低 。 
关键 词 : 冰川 反照 率 ; 黑河 流域 ，MOD10A1; Landsat TM/ETM+; 时 空 分 布 


冰川 反照 率 作为 影响 辐射 收文 最 为 敏感 的 要 
素 ,控制 着 大 范围 内 融 水 产 出 的 时 空 分 布 ,是 冰川 
消融 模拟 的 关键 参数 ,其 波动 变化 将 直接 影响 冰川 
能 量 一 物质 平衡 的 变化 "。 冰 川 反照 率 对 于 气候 变 
化 极其 敏感 ,在 不 同 的 时 间 与 空间 上 ,冰川 反照 率 
的 差异 较 大 了。 在 全 球 气候 变 暧 的 背景 下 ,监测 和 
计算 冰川 反照 率 对 于 估算 冰川 消融 意义 重大 AG 
统 的 冰川 反照 率 人 研究 多 以 气象 站 和 地 物 波谱 仪 为 
基础 ,具有 时 空 尺度 受 限 观测 难度 大 等 局 限 性 呈 。 
随 着 遥感 技术 的 不 断 发 展 ,获得 高 时 空 分 辨 率 .长 
时 间 序 列 的 冰川 反照 率 数据 成 为 了 可 能 ,为 冰川 反 
照 率 研究 提供 了 新 的 技术 手段 1。 

MOD10A1 逐日 反照 率 产 品 具有 较 高 的 时 空 分 
辨 率 , 国 内 外 学 者 对 该 套数 据 进 行 了 较 多 的 应 用 和 
研究 。 已 有 研究 表明 ,在 格陵兰 冰 盖 地 区 MODIS 逐 
日 反照 率 产 品 可 以 反映 一 般 的 季节 性 变化 , 且 Terra 
卫星 的 反 演 结果 略 优 于 Aqua B.E cde ES RD S 
在 土耳其 山地 冰川 区 ,MODIS 逐 日 反照 率 产 品 与 实 
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地 测量 结果 较为 一 致 (大 约 高 估 了 10% ) ,同时 ,在 
高 海拔 处 的 适用 性 更 她。 基于 MODI10A1 和 逐日 反 
照 率 产品 的 较 好 适用 性 ,可 依据 该 套数 据 分 析 单 条 
冰川 .甚至 大 范围 冰川 区 的 反照 率 时 空 分 布 特征 及 
变化 规律 号, 还 可 以 利用 MODI10A1 产品 为 驱动 ， 
重建 冰川 年 兆 物 质 平衡 ”。 前 人 对 于 MOD10A1 数 
据 的 应 用 大 多 着 眼 于 大 尺度 的 区 域 或 针对 小 尺度 
某 一 冰川 进行 精度 验证 和 时 空 变化 研究 ,而 对 于 中 
尺度 流域 的 冰川 反照 率 关 注 较 少 。 本 文 以 黑河 上 
游 典 型 冰川 (七 一 冰川 . 羊 龙 河 五 号 冰川 . 八 一 冰 
川 .二 一 冰川 ) 为 研究 对 象 ,利用 MOD10A1 产品 研 
究 2000 一 2018 年 反照 率 时 空 变化 特征 ,以 期 为 定量 
计算 冰川 消融 提供 依据 。 


1 研究 区 概况 


黑河 流域 位 于 祁连山 中 段 北 坡 (38° ~ 42°N, 
98? ~ 101°30'E) ,是 我 国 第 二 大 内 陆 河 ”"。 根 据 
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第 二 次 冰川 编目 资料 ,黑河 流域 共 发 育 有 冰川 952 
条 ,总 面积 达 293.60 km" "^, 

本 文选 取 黑 河流 域 上 游 的 4 条 典型 冰川 (图 1)， 
自 西向 东 分 别 为 托 来 山北 坡 的 七 一 冰川 , 素 珠 链 峰 
北 坡 的 羊 龙 河 五 号 冰川 , 玻 款 山 南 坡 的 八 一 冰川 和 
萌 芦 沟 流 域 源 头 的 十 一 冰川 时。4 条 典型 冰川 的 
面积 分 别 为 2.53 km?、5.68 km?,2.39 km’ Fil 0.50 km’; 
海拔 高 度 分 别 在 4323.70 ~ 5114.70 m、4291.10 ~ 
5406.70 m ,4567.50-4823.70 m 和 4338.40 ~ 4774.60 
m [X [B] j^ 
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图 1 研究 区 地 理 位 置 图 


Fig. 1 The location of the four representative glaciers 


in Heihe River Basin 


2 数据 及 方法 


2.1 数据 源 

2.1.1 MOD10A1 反 照 率 产品 MOD10A1 是 美国 航 
空 航 天 局 (NASA) 陆 地 产品 组 研制 的 三 级 逐日 积 雪 
数据 产品 ,空间 分 辩 率 为 500 m。 本 文 从 美国 雪 冰 
数据 中 心 (NSIDC ,https:/nsidc.org/) 获 取 了 产品 编码 
y h25v05 , JA, H 8] 2g 2000 4E. 2. H $1 2018 Æ 12 H 
的 影像 数据 。 

2.1.2 Landsat TM/ETM + iÉ E X5 4% Landsat TM/ 
ETM+ 系 列 遥 感 影像 来 自 于 美国 地 质 调查 局 影像 数 
据 库 中 心 (USGS https://earthexplorer.usgs.gov/) 。 选 
取 与 下 文 介绍 气象 站 实测 数据 获取 日 期 对 应 的 、 云 
和 积 雪 较 少 的 6 景 TMETM+ 影 像 ,具体 影像 信息 如 
表 1 所 示 。 

2.1.3 反照 率 实测 资料 七 一 冰川 自动 气象 台 站 实 
测 反 照 率 数据 来 自前 人 研究 ,获取 了 2006 年 8 
月 1 日 .8 月 25 日 .9 月 18 日 ,2007 年 8 月 12 日 .9 月 


ChinaXiv 合 作 期 刊 


6j 王 俊 瑶 等 :基于 MOD10AI 的 祁连山 黑河 流域 典型 冰川 反照 率 时 空 变化 研究 1397 


表 1 Landsat TM/ETM+ 影 像 信息 
Tab.1 Information of Landsat TM/ETM- images 


获取 时 间 行列 号 fede THIER 
2006-08-01 135/033 TM 30 
2006-09-18 135/033 TM 30 
2007-09-21 135/033 TM 30 
2006-08-25 135/033 ETM+ 30 
2007-08-12 135/033 ETM+ 30 
2007-09-13 135/033 ETM+ 30 


13 日 .9 月 21 日 共 6 日 的 数据 作为 验证 数据 ,七 一 冰 
JI A RAH (9745.2 E, 39° 14.7 N ,海拔 4473 m, 
坡度 10") 配 设 有 CNR1 净 辐 射 表 ,安装 高 度 为 1.5 
m, 上 且 与 水 平面 保持 平行 。 

十 一 冰川 反照 率 实测 数据 来 自 于 HandHeld™an 
地 物 波谱 仪 实际 测量 ,仪器 观测 波长 精度 为 +/-1 
nm@700 nm。 本 研究 所 用 数据 测量 时 间 为 2013 年 8 
月 13 日 一 8 月 31 日 每 天 的 上 午 9:00 至 下 午 13:00， 
天 气 状况 为 晴空 .光照 条 件 稳 定 ;测量 地 点 为 十 一 
冰川 消融 区 花 杆 点 处 ,测量 时 探头 与 冰 面 距离 为 
1.3 m, 测 量 时 每 10 ~ 15 min 对 参考 白板 进行 优化 ， 
利用 标准 太阳 入 射 辐射 光谱 计算 的 方法 将 ASD 观 
测 的 反射 光谱 转换 为 反照 率 。 由 于 对 应 时 间 段 内 
MOD10A1 产 品 反照 率 有 效 值 缺 失 , 所 以 该 数据 用 于 
与 测量 日 期 前 后 的 MOD10A1 产品 反照 率 进行 对 比 
分 析 。 
2.1.4 气象 再 分 析 资 料 ERA-Interim 再 分 析 资 料 
来 自 欧洲 中 期 数值 预报 中 心 (ECMWE https://www. 
ecmwf.int/) ,本 文选 取 2000—2018 年 0.75? 4) SES AY 
月 平均 气温 、 总 降水 和 降雪 数据 。HAR (The High 
Asia Refined analysis) 来 自 柏林 工业 大 学 (https:// 
www.klima.tu-berlin.de/) ,本 文选 取 2002 年 云 量 因子 
日 均值 数据 。 
2.1.5 DEM 数字 高 程 模型 (DEM) 是 美国 航空 航天 
Jay (NASA ) 和 美国 国防 部 国家 测绘 局 (NIMA ) 等 多 部 
门 共 同 测量 的 SRTM v4.1, 该 数据 获取 自 地 理 空间 
数据 云 (http://www.gscloud.cn/) ,用 来 辅助 Landsat 
TM/ETM+ 数 据 的 地 形 校正 分 析 。 
2.2 方法 

本 文 利 用 实测 数据 进行 Landsat TM/ETM+ 反 演 
反照 率 和 MOD10A1 产品 反照 率 的 误差 分 析 ,批量 
提取 MOD10A1 产 品 在 冰川 区 内 有 效 值 进行 反照 率 
日 均值 的 计算 ,并 分 析 2000 一 2018 年 黑河 流域 典型 
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冰川 反照 率 的 时 空 分 布 特征 9 技术 路 线 如 图 2 所 示 [e] a= Lok ~ Ln (2) 

2.2.1 MODIOALT& RE JE > ese! = MODIOALT Jc FR 255 
B7 Lus 3) 


率 数据 层 中 积 雪 反照 率 的 数值 介 于 0 ~ 100( #22). 
MOD10A1 原始 产品 进行 格式 及 投影 转换 ,并 依据 第 
二 次 冰川 编目 数据 掩 膜 提取 冰川 区 MOD10A1 产品 
像 元 值 后 ,进行 有 效 值 提 取 和 计算 。 
| MOD10A1 产 品 | 
( RAE ) 


Y 
按 冰 川 区 掩 膜 


有 效 值 提取 


Y 
MODIS 逐 日 反照 率 


Y 
误差 分 析 
| 


气象 站 实测 反照 率 | 


Landsat TM/ETM+ 
反 演 反照 


Y 
黑河 流域 典型 冰川 


黑河 流域 典型 冰川 i 
反照 率 时 间 变 化 


反照 率 空间 分 布 


图 2 技术 路 线 


Fig. 2 Technical chart 


表 2 MOD10A1 产 品 中 像 元 值 的 含义 
Tab. 2 Meaning of the pixel value of MOD10A1 


像 元 值 含义 像 元 值 含义 
0-100 ” 积 雪 反 照 率 150 云层 覆盖 

101 不 确定 的 地 物 类 型 | 251 自身 阴影 

111 夜间 252 陆地 界限 不 匹配 
125 无 积 雪 覆盖 的 陆地 | 253 BRDF 无 效 

137 内 流 湖 254 没有 匹配 的 数据 
139 开阔 水 体 250 缺失 的 数值 


2.2.2 Landsat TM/ETM+ 反 照 率 反 演 算法 
(1) 辐射 定 标 

按照 遥感 影像 头 文件 说 明 ,将 影像 的 原始 DN 

值 转换 为 具有 物理 意义 的 大 气 顶 部 辐射 亮度 值 L 

(Weems sr!) ,辐射 定 标 公式 为 : 

L=a:DN+B 


(1) 


式 中 :a.B6 分 别 表示 辐射 定 标 系数 中 增益 、 偏 移 量 ， 
Li 和 Li 分 别 为 辐射 亮度 值 的 动态 变化 范围 。 

(2) 地 形 校正 

山区 由 于 地 形 复杂 ,受到 遮蔽 作用 的 影响 , 阴 
了 昌 坡 冰雪 地 物 在 影像 中 的 辐射 亮度 产生 明显 差 
异 。 因 此 必须 通过 地 形 校正 减少 山区 地 形 遮 蔽 作 
用 对 影像 辐射 亮度 的 影响 ,提高 反 演 地 表 反 射 率 精 
RE, AS x fii Hd ENVI Topographic Correction 工具 
选择 SCS+C 模型 进行 地 形 校正 ,此 模型 引入 半 经 验 
的 C 系 数 来 说 明 散 射 辐射 ,对 SCS 模 型 进行 改进 ”， 
具有 调节 过 度 校 正 的 作用 2 。 

(3) 大 气 校正 

电磁 波 在 穿 过 大 气 层 时 会 受到 大 气 衰 弱 的 影 
响 ,造成 遥感 反 演 反射 率 的 误差 ”。 需 要 使 用 大 气 
校正 来 消除 影响 ,将 大 气 层 顶 部 的 辐射 量 转换 为 反 
映 地 物 真 实 信息 的 地 表 反 射 率 。 本 文选 用 FLAASH 
模型 进行 大 气 校正 ,输入 参数 如 表 3 所 示 。 

(4) 各 向 异性 校正 

由 于 地 表 反 照 率 是 各 观测 角 的 反射 率 的 积分 ， 
而 卫星 仅 在 某 一 个 或 几 个 角度 进行 观测 ,因而 必须 
应 用 BRDF 模 型 进行 各 向 异性 校正 以 减 小 反射 率 差 
异 ”。 本 人 研究 采用 C 因子 法 BRDF 校正 ,此 方法 利 
用 全 球 MODIS BRDF 产品 中 约 155 亿 像 元 值 拟 合 得 
到 的 BRDF 模 型 进行 校正 ”。 

(5) 罕 波 段 反射 率 转化 为 宽 波 段 反照 率 

由 于 卫星 迁 感 测量 的 是 罕 波 段 不 连续 的 反射 
率 , 必 须 将 各 个 鹤 波 段 的 光谱 反射 率 转化 为 宽 波段 
反照 率 。 本 文采 用 Liang 等 汪 提 出 的 罕 波 段 转 宽 波 
段 的 转化 公式 : 

a=0.443a, + 0.317a, + 0.240a, (4) 

式 中 :a 表示 反照 率 反 演 值 ,a,a;,a; 分 别 表示 第 1、 
2.3 波 段 反 射 率 。 


表 3 FLAASH 大 气 校正 输入 参数 (以 2006082SETM+ 数 据 为 例 ) 
Tab. 3 Parameters required input into FLAASH for atmospheric correction (taken 20060825 ETM+ data as an example) 


影像 中 心 点 传感器 类 型 飞行 日 期 海拔 高 度 /km 像 元 大 小 大 气 模式 气 溶胶 模式 
38°54'43.79"N 
_08- i x < 
ost ae SUPE Landsat ETM+ 2006-08-25 4.773 30 mx30 m MLS 乡村 
初始 能 见 度 /km RI scale fey EE/km CO0: 混 合 比例 /ppm MODTRAN 分辨 率 /cm 散射 模型 DISCORT 流 量 
33 1.25 370 15 scale DISCORT 8 
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2.2.3 误差 分 析 黑河 流域 的 典型 冰川 受 复杂 地 形 
和 风 吹 雪 的 影响 ,空间 异 质 性 较 大 ”“。 自 动 气 象 站 
测定 的 有 效 视野 为 10 mx10 m 左 右 ,而 MOD10A1 反 
照 率 产品 的 空间 分 辨 率 为 500 mx500 m, 两 者 观测 视 
野 间 的 差距 较 大 ,必然 会 导致 冰川 反照 率 的 差异 。 
因而 ,在 进行 误差 分 析 时 ,必须 要 考虑 空间 尺度 效 
应 对 结果 的 影响 。 本 文采 用 Liang 提 出 的 利用 高 空 
间 分 状 率 数据 产品 ,将 地 面 实测 数据 升 尺度 到 MO- 
DIS 空间 尺度 进行 误差 分 析 。 

精度 验证 主要 根据 实测 的 反照 率 将 Landsat 
TMETM+ 反 演 反 照 率 瞬时 值 转换 为 日 值 , 并 对 转换 
后 的 Landsat TWETM+ 影 像 建立 17x17 像 元 的 滑动 
平均 窗口 ,使 其 空间 分 辨 率 上 升 到 510 mx510 m, £f 
对 其 重 采样 至 500 mx500 m, 并 与 MOD10A1 反 照 率 
关 品 进行 空间 配 准 。 通 过 得 到 的 Landsat TM/ETM+ 
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反 演 反照 率 结果 对 MODI0A1 反 照 率 产品 进行 逐 像 
元 精度 验证 。 评 价 指标 包括 相关 系数 (Correlation 
Coefficient, r) 、 绝 对 误差 (Absolute Error, AE) .平均 
绝对 误差 (Mean Absolute Error, MAE ) 和 均 方 根 误差 
(Root Mean Square Error, RMSE )。 


3 结果 与 分 析 


3.1 MODIS 反照 率 产品 在 黑河 流域 典型 冰川 的 误 
差分 析 

对 MOD10A1 反 照 率 产 品 提取 结果 与 七 一 冰川 
表面 自动 气象 站 测 得 的 实测 反照 率直 接 进行 对 比 
(Fl 3). MODIOAI 反照 率 产品 提取 结果 与 七 一 冰 
川 表面 自动 气象 站 测 得 的 实测 反照 率 绝 对 误差 (4E) 
在 -0.26~0.08 之 间 , 平 均 绝 对 误差 (MAE) 为 -0.09， 
均 方 根 误差 (RMSE ) 为 0.15 ,两 者 随时 间 的 变化 趋势 


1.07 (p ~ MODI0A1 产 品 反照 率 

一 一 实测 反照 率 

0.8 上 — MOD10A1 产 品 反照 率 线性 趋势 
— 实测 反照 率 线性 趋势 


图 3 MODI10A1 产 品 反 照 率 与 实测 反照 率 比较 
Fig. 3 Comparison of MOD10A1 glacier albedo and ground measured albedo 


较为 一 致 。BRDF 模型 的 误差 IRA ROCA a Tu 
等 原因 是 MOD10A1 产品 出 现 误 差 的 主要 原因 ” 。 

进而 利用 Liang 提出 的 方案 ,用 Landsat TM/ 
ETM+ 影 像 反 演 反照 率 对 MOD10A1 反 照 率 产品 进 
行 误差 分 析 。 根 据 气象 站 的 日 反照 率 数据 与 其 对 
应 日 期 的 Landsat TMETM+ 影 像 反 演 出 的 瞬时 反照 
率 数 据 进 行 线性 回归 分 析 , 计算 七 一 冰川 单条 冰川 
的 日 反照 率 值 (图 4)。 

以 MOD10A1 反 照 率 产品 在 七 一 冰川 的 适用 性 
来 代表 该 产品 在 黑河 流域 4 条 典型 冰川 的 适用 性 。 
Landsat 反 演 反 照 率 与 MOD10A1 产 品 反照 率 的 对 比 结 
果 表 明 , 两 者 绝对 误差 (AE) 在 -0.06 ~ 0.09 之 间 ,平均 
绝对 误差 (MAE) 为 0.01, 均 方 根 误差 (RMSE) 为 0.05。 
表明 该 产品 在 黑河 流域 具有 较 好 的 适用 性 (图 5)。 


y=0.94x+0.01 
3 


0.1 0.3 0.5 0.7 
Landsat TM/ETM-- 影像 反 演 反照 率 瞬 时 值 


图 4 Landsat TM/ETMT+ 反 演 反照 率 瞬 时 值 与 冰川 表面 实 
测 反照 率 日 值 的 散 点 关系 
Fig. 4 The scatter relationship between Landsat TM/ETM+ 


derived albedo and ground measured albedo 
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(a) Landsat TM/ETM-- 反 演 反照 率 (b) MOD10A1 产 品 反 演 反照 率 
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Z 图 例 e Z 图 例 Z 
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umm 0.21~0.25 m 0.21~0.25 
mu 0.26~0.30 mu 0.26~0.30 
= 0.31~0.35 = 0.31~0.35 
0 1 km 0 1 km 
97°45'E 97°46E 97°45'E 97°46'E 
O57 (c) 散 点 关系 
* 
E 
= 0.4 上 
LN 
E 
A 03 
S 
E 
0.2 ， 
0.2 0.3 0.4 0.5 
Landsat TM/ETM-- 影像 反 演 反照 率 
图 5 Landsat TMETM+ 影 像 反 照 率 空间 尺度 转换 结果 与 MOD10A1 产 品 反照 率 
Fig. 5 Albedo derived from Landsat TM/ETM+ images and MOD10A1 albedo product 
3.2 黑河 流域 典型 冰川 反照 率 时 间 变 化 特征 di 
分 析 MODIOA1 产品 反照 率 的 年 内 变化 结果 | 
(图 6), 八 一 冰川 反照 率 在 3 月 有 小 幅 上 升 ,5 月 后 
迅速 下 降 在 8 月 达到 最 低 值 ,随后 大 幅度 增 大 并 于 i 0.5. 
11 月 到 达 最 高 值 ,12 H RAE 3. H CHR 5 04085) e i 
少 的 趋势。4 一 5 月 可 能 由 于 降雪 频率 有 所 增加 , 冰 x 04r 
川 表面 反照 率 有 小 幅度 上 升 ;6 一 8 月 ,冰雪 消融 作 T 
用 最 强 , 反 照 率 在 此 时 最 低 ;9 一 11 月 , 随 着 冰川 消 l 
融 作 用 减弱 ,反照 率 不 断 升 高 ,并 达到 峰值 。 0.2 
0 2 4 6 8 10 12 


黑河 流域 4 条 典型 冰川 反照 率 随 时 间 变 化 趋势 
较为 一 致 (图 7) ,最 低 值 均 出 现在 2013 年 ;七 一 冰 
川 . 羊 龙 河 五 号 冰川 和 八 一 冰川 最 大 值 均 出 现在 
2008 年 ,虽然 十 一 冰川 反照 率 最 大 值 出 现在 2003 
年 ,但 2008 年 冰川 反照 率 是 次 高 峰 , 日 与 最 大 值 的 差 
异 较 小 。4 条 典型 冰川 反照 率 在 2000 一 2018 年 均 呈 
现 微弱 的 下 降 趋 势 , 年 均 下 降 速 率 分 别 为 -0.0004、 
-0.002 -0.0004 和 -0.0008 ,其 中 , 羊 龙 河 五 号 冰川 
反照 率 下 降 速度 最 快 ,十 一 冰川 反照 率 下 降 速率 次 
之 。4 条 冰川 均 未 通过 95% 的 置信 度 检 验 ,说 明 4 条 


6 2000 一 2018 年 八 一 冰川 多 年 平均 反照 率 逐 月 变化 
Fig.6 Monthly variations of albedo on the Bayi Glacier 
from 2000 to 2018 


冰川 年 平均 反照 率 的 下 降 趋势 不 显著 。 而 根据 


ERA- Interim 提取 的 总 降水 和 降雪 数据 ,年 均 固 态 


降水 与 液态 降水 比例 下 降 速率 分 别 为 -0.01、 


-0.02、-0.01 和 -0.02。 羊 龙 河 五 号 冰川 和 十 一 冰川 


年 均 降 水 固 液 比 下 降 速 度 最 快 ,这 可 能 是 导致 这 两 
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图 7 黑河 流域 4 条 典型 冰川 反照 率 的 年 际 变化 (MOD10A1) 
Fig.7 Annual variations of albedo on the four representative glaciers in Heihe River Basin (MOD10A1) 


个 冰川 反照 率 下 降 速率 最 快 的 原因 。 川 的 反照 率 空 间 分 布 为 例 ( 图 8) ,两 条 冰川 均 表现 
3.3 黑河 流域 典型 冰川 反照 率 空间 变化 特征 出 冰川 区 下 部 反照 率 低 于 上 部 反照 率 的 特征 ,但 


为 了 反映 黑河 流域 典型 冰川 反照 率 的 空间 变 。 八 一 冰川 下 部 反照 率 与 上 部 反照 率 空间 差异 
化 特征 ,本 文 依据 二 次 冰川 编目 的 冰川 中 值 面积 高 。 较 小 。 


度 , 提 取 对 应 高 度 线 作 为 划分 冰川 区 上 部 和 下 部 的 已 有 研究 表明 ,每 年 7 一 8 月 为 黑河 流域 典型 冰 
依据 。 川 的 消融 旺季 ”3,4 条 冰川 多 年 反照 率 均 值 及 消融 
以 2008 年 8 月 14 日 羊 龙 河 五 号 冰川 \ 八 一 冰 旺季 反照 率 均 值 如 图 9 所 示 。 年 平均 反照 率 八 一 冰 
98°33'E 98?34'E 98?35'E 98°53'E 98?54'E 
(a) 羊 龙 河 五 号 冰川 z| 四 从 一 冰川 " 
=| $ 
zx z &| 冰川 区 下 部 一 A 
>t E 
$ S 
z z 
z z al JS 
x 图 例 IE 3 
© | \ 羊 龙 河 五 号 冰川 a 5 图 例 RR a 
- 100 m 等 高 线 八 一 冰川 N: 
E TU. : 冰川 上 下 部 分 界线 us m 等 高 线 一 冰川 上 下 部 分 界线 
Z0 1 km ma 0.31-0.45 mm 0.61~0.75 | Z mm 0.31~0.45 mm 0.61~0.75 Q 0.5km 
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图 8 羊 龙 河 五 号 冰川 和 八 一 冰川 反照 率 空间 差异 (MOD10A1) 
Fig. 8 The spatial difference of albedo on Yanglonghe glacier No. 5 and Bayi Glacier (MOD10A1) 
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图 9 黑河 流域 4 条 典型 冰川 反照 率 值 空间 分 布 

(MODI0A1) 


Fig.9 Spatial distribution of albedo on the four representative 
glaciers in Heihe River Basin (MOD10A1) 


川 最 高 . 羊 龙 河 五 号 冰川 次 之 .七 一 冰川 次 低 .于 一 
冰川 最 低 ; 消融 旺季 冰川 反照 率 羊 龙 河 五 号 冰川 最 
高 .七 一 冰川 次 之 、 八 一 冰川 次 低 、 十 一 冰川 最 低 ， 
消融 旺季 4 条 冰川 反照 率 基 本 呈现 随 海拔 上 升 而 上 
升 的 特征 ,但 整个 流域 是 否 存在 同样 的 冰川 反照 率 
随 着 海拔 高 度 而 变化 的 规律 需 在 日 后 的 研究 中 加 
以 探讨 。 


4 讨论 


4.1 云 量 对 于 冰川 反照 率 的 影响 

大 气 中 的 云 量 是 影响 反照 率 的 重要 因素 ” ,本 
文 判 定 云 量 小 于 0.3 时 为 晴天 ,大 于 0.8 时 为 阴 天 ， 
介 于 两 者 之 间 为 多 云 ””"。2002 年 不 同 云 量 因 子 下 4 
条 冰川 反照 率 变化 如 表 4 所 示 , 阴 天 时 冰川 反照 率 
会 比 晴 天 时 有 所 增加 。 主 要 原因 是 云 对 于 近 红 外 
波段 的 辐射 吸收 较 强 ,了 明天 云 吸收 了 较 多 的 近 红外 
波段 辐射 ,到 达 冰 川 表面 的 可 见 光 波段 占 整 个 太阳 
入 射 短波 中 的 比例 增 大 。 同 时 冰川 表面 对 于 可 见 
光波 段 辐射 的 反照 率 最 高 ,因而 阴 天 下 冰川 反照 率 


表 4 不 同 云 量 因子 下 反照 率 变化 
Tab.4 Albedo variations under different 


cloud cover factors 


冰川 名 称 晴天 下 冰川 反照 率 ” 阴 天 下 冰川 反照 率 
七 一 冰川 0.35 0.36 
羊 龙 河 五 号 冰川 0.39 0.50 
八 一 冰川 0.45 0.48 
十 一 冰川 0.31 0.41 
平均 值 0.38 0.44 


增加 ,而且 辐射 在 云 底 和 冰川 表面 之 间 多 次 反射 也 
会 加 剧 冰川 反照 率 的 增加 效应 。 
4.2 气候 因子 对 冰川 反照 率 的 影响 

降水 类 型 和 气温 变化 是 影响 冰川 反照 率 变化 
的 重要 影响 因素 ”。 根 据 ERA-Interim 再 分 析 资 料 
计算 的 数据 显示 ,黑河 流域 4 条 典型 冰川 2000 一 
2018 年 4 条 冰川 年 平均 气温 呈 微 弱 上 升 趋势 年 平 
均 固 液 比 均 呈 微弱 的 下 降 趋 势 , 气 温 趋 势 与 反照 率 
变化 趋势 呈 负 相关 (相关 系数 + 分 别 为 -0.37、 
-0.17、-0.38 和 -0.66) 、 固 液 比 趋势 与 反照 率 变 化 趋 
势 呈 正 相 关 ( 相 关系 数 > 分 别 为 0.22.0.19 .0.23 和 
0.12)。 不同 时间 段 内 ,消融 季 气 温和 固 液 比 对 于 反 
照 率 的 影响 程度 不 尽 相 同 ,可 能 与 其 他 原因 相关 ， 
如 冰川 反照 率 与 积 雪 密度 和 积 雪 污 化 物 颗粒 浓度 
大 致 呈 反 相关 的 线性 关系 3”, 年 平均 气温 、 年 总 降 
水 量 、 黑 碳 、 粉 侍 沉 积 状况 等 也 不 同 程度 的 影响 着 
冰川 反照 率 。 


5 结论 

通过 对 Landsat TMETM+ 影 像 和 2000 一 2018 年 
MOD10A1 产 品 提供 的 反照 率 数 据 进行 分 析 , 研 究 了 
黑河 流域 4 条 典型 冰川 反照 率 的 时 空 分 布 和 变化 特 
征 ,得 出 以 下 结论 : 

(1) MOD10A1 反 照 率 产 品 提取 结果 与 七 一 冰 
川 表面 自动 气象 站 实测 的 反照 率 随 时 间 的 变化 趋 
势 一 致 , 绝 对 误差 (AE) 在 -0.26 ~ 0.08 之 间 ,平均 绝 
对 误差 (MAE) 为 -0.09, 均 方 根 误差 (RMSE ) 为 
0.15。MOD10A1 反 照 率 产品 与 Landsat TM/ETM+ 影 
像 反 演 反 照 率 值 差异 较 小 ,绝对 误差 (AE) 在 -0.06 ~ 
0.09 之 间 ,平均 绝对 误差 (MAE ) 为 0.01 , 均 方 根 误差 
(RMSE ) 为 0.05 ,证 明 MODI10OA1I 反 照 率 产 品 在 黑河 
流域 具有 较 好 的 适用 性 。 

(2) 黑河 流域 4 条 典型 冰川 反照 率 随 时 间 变 化 
较为 一 致 ,最 低 值 均 出 现在 2013 年 ;七 一 冰川 . 羊 龙 
河 五 号 冰川 和 八 一 冰川 最 大 值 均 出 现在 2008 年 , 虽 
然 十 一 冰川 反照 率 最 大 值 出 现在 2003 年 ,但 2008 
年 冰川 反照 率 是 次 高 峰 ,日 与 最 大 值 的 差异 较 小 。 
4 条 典型 冰川 反照 率 在 2000 一 2018 年 均 呈 微弱 的 下 
降 趋势 。 

(3) 羊 龙 河 五 号 冰川 和 八 一 冰川 2 条 冰川 均 呈 
冰川 区 上 部 冰川 反照 率 高 于 下 部 反照 率 的 特征 , 八 
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一 冰川 上 下 部 反照 率 空间 差异 较 小 。 1052. 
(4) 年 平均 冰川 反照 率 八 一 冰川 最 高 , 羊 龙 河 [10] $5734, 蒙 吉 军 , 李 枫 . 基于 空间 异 质 性 的 黑河 中 游 水 资源 脆弱 
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Spatiotemporal variation of albedo of four representative glaciers in the Heihe 
River Basin based on multi-source data 


WANG Jun-yao, HUAI Bao-juan, WANG Ye-tang SUN Wei-jun, ZHANG Wu-ying 
(College of Geography and Environment, Shandong Normal University, Jinan 250300, Shandong, China) 


Abstract: Glacier albedo is a critical parameter that controls the process of radiation budget, the fluctuation of 
which directly affects the change in glacier energy—mass balance. By remote sensing, studies on the characteristics 
of glacier albedo and its spatiotemporal distribution can be taken as an important reference for estimating glacier 
melting and downstream river runoff. In this study, albedo dataset retrieved from Landsat TM/ETM+ images and 
MOD10A1 daily albedo product were used as the main data sources, and Landsat TM/ETM+ images and MOD10A1 
data’ error were analyzed by field measurements. Then, the spatiotemporal variations of albedo were investigated on 
four representative glaciers in Heihe River Basin of MOD10A1 from 2000 to 2018. The results showed that there 
was a certain bias between albedo value derived from MODIOA1 products and ground measured, but the time series 
was similar; there is a difference in inversion results between MOD10A1 and Landsat TM/ETM+, suggesting a good 
applicability of glacier albedo in Heihe River Basin. On the interannual scale, the albedo of four glaciers was 
decreasing. Among the spatial scale, the annual average albedo of Bayi glacier was the highest, followed by 
Yanglonghe glacier No. 5, Qiyi glacier, and Shiyi glacier; the average albedo of Yanglonghe glacier No. 5 was the 
highest in the ablation season, followed by the Qiyi glacier, Bayi glacier, and Shiyi glacier. 


Keywords: glacier albedo; Heihe River Basin; MOD10A1; Landsat TM/ETM+; temporal and spatial variation 


